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ABSTRACT 

The fraction of a fireball kinetic kinetic energy that can be radiated away by internal shocks is sensitive to the 
amplitude of initial fluctuations in the fireball. We give a simple analytical description for dissipation of modest- 
amplitude fluctuations and confirm it with direct numerical simulations. At high amplitudes, the dissipation occurs 
in a non-linear regime with efficiency approaching 100 %. Most of the explosion energy is then radiated away by 
the prompt GRB and only a fraction remains to be radiated by the afterglow. 

Subject headings: gamma-rays: bursts 



1. INTRODUCTION 

Cosmological gamma-ray bursts (GRBs) are huge explosions 
of energy ~ 10 53 ergs, which may be triggered, e.g., by coales- 
cence of neutron stars and/or black holes. The created fireball 
expands relativistically, with a Lorentz factor Y ~ 10 2 . In the 
most likely scenario, the observed gamma-rays are produced at 
times t ~ 10 2 — 10 5 s after the explosion, when the relativistic 
outflow gets optically thin and before it is decelerated by the 
surrounding medium (see, e.g., Piran 1999 for a review). The 
outflow has a fluctuating velocity profile, and the 7-rays are 
generated by internal shocks that develop when faster shells try 
to overtake slower ones (Rees & Meszaros 1994). 

The internal dissipation was discussed and simulated numeri- 
cally in a number of works (e.g., Kobayashi, Sari, & Piran 1997; 
Daigne & Mochkovitch 1998; Lazzati, Ghiselhni, & Celotti 
1999; Panaitescu, Spada, & Meszaros 1999; Spada, Panaitescu, 
Meszaros 1999; Kumar 1999). It was concluded that only a few 
percent of the total kinetic energy is emitted to infinity (e.g., 
Panaitescu et al. 1999; Kumar 1999). Then the energy of the 
prompt GRB should be dominated by the early afterglow asso- 
ciated with an external shock. Observations show the opposite. 
Besides, the low efficiency requires huge explosion energies. 

The conclusion that the internal dissipation has a low effi- 
ciency was based on certain assumptions: 

(1) A specific probability distribution was assumed for the out- 
flow Lorentz factor, such that the root-mean-square (rms) of the 
fluctuations was less than 1 / y3 (see eq. [11]). In fact, the rms 
can be much higher. 

(2) The fluctuations were taken with a white (flat) spectrum. 

(3) It was often assumed that the dissipated energy is equally 
distributed between the electrons, protons, and magnetic field, 
and only the electron part is radiated. In fact, the radiative ca- 
pability of shocked matter is highly uncertain and dependent on 
poorly understood plasma processes. The high radiative capa- 
bility is favored by the observed high luminosities of the prompt 
GRBs. If most of the dissipated energy is transferred to the 
electrons and radiated, the outflow stays at low pressure and the 
coasting matter gets concentrated into thin shells (caustics). 
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In § 2, we give a simple analytical description for dissipation 
of white fluctuations with modest rms (the linear regime). Then 
the efficiency is proportional to the mean square of the initial 
fluctuations. In §3, we study high-amplitude fluctuations. The 
efficiency then approaches 100 %. The analytical results are 
illustrated with numerical simulations. 

2. LOW-AMPLITUDE FLUCTUATIONS 

The size of the central engine, r- m ~ 10 7 — 10 8 cm corre- 
sponds to a typical time scale t ln = r- ln /c ~ 1 ms. The outflow 
generated by the engine can have highly correlated parameters 
on this time scale, so that portions At i m can be considered 
as individual shells in the continuous outflow. The outflow is 
then discretized as a sequence of such shells, i = 1, N, with 
mass m, and energy input e,-. One after another, the shells are 
accelerated to a Lorentz factor r,- = e,/ rrii (e, gets converted 
into the bulk kinetic energy) on the time-scale Yiti n < 1 s (see 
Piran 1999). Subsequent evolution of the outflow proceeds in 
the coasting regime until the shells begin to collide. 

Suppose Yi is a white noise with the average value T av and 
initial rms r rms0 < r av . Masses of the shells can be taken 
equal, m$ = mo, or also fluctuating - the results will be the 
same. The initial scale of the fluctuations is Ao ~ Hn- The 
collision process starts at 



to 



Ao 



r„ 



(i) 



After time t ~ 2to, the first generation of shell coalescence 
has been done (the typical mass of a shell increases by a fac- 
tor of two), then the second generation occurs, and so on. In 
the process of hierarchical coalescence, the scale of the fluctu- 
ations increases, A = A(i), as well as the average mass of an 
individual shell, X(t)/Xo — m(t)/mo = K(t). 

The fluctuations look especially simple when viewed from 
the frame moving with Lorentz factor T av . We will denote 
quantities measured in this frame by symbols with tilde. The 
fluctuation velocity is v,-/c = (Yf - r 2 v )/(r 2 + T 2 V ). Given 
r rms < r av , we have fyfc « (I\ - r av )/T av with the average 
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and the rms, 



(2) 



The outflow motion is thus decomposed into two parts: the rela- 
tivistic motion of the center-of-momentum (CM) with a Lorentz 
factor Tcm ~ Tav and superimposed non-relativistic fluctua- 
tions. 

The hierarchical coalescence of shells can be expected to oc- 
cur in a self-similar regime, so that the collisional time-scale 
for one generation is about the time passed since the explosion, 



A 



t = 



t 



r c 



M 



(3) 



Here, A = Tom A is the scale of the fluctuations in the CM- 
frame. After coalescence of K shells with initial momentum 
po <~ moSrmso, the new big shell has a momentum p <~ \f~Kp~o 
(this is the random walk formula: K momenta ~ po are 
summed with random sign). We thus have, 

rm) rms ~ VKm v rms0 = J— m v rms0 . (4) 
V m o 

Combining (2), (3) and (4), we get the self-similar solution de- 
scribing the hierarchical coalescence at t > to, 



to 
mo 



A_ 



to J 



2/3 



rmsO 



-1/3 



(5) 



Note that the collisions establish a Gaussian distribution of v. 
One can therefore prescribe a temperature to the system of 
shells, T, which is related to the average kinetic energy by 
mv 2 ms /2 = kT/2 (shells have one degree of freedom). From 
equation (5) it follows that T = const. The evolution of shells 
can thus be described as isothermal sticking together. 
The free energy of the outflow is given by 



volume increases proportionally to t 2 . As a result of adiabatic 
cooling, the volume-integrated radiation energy is reduced as 
£-2/3 jjjg ener gy conservation law implies that the adiabatic 
cooling is accompanied by a regular radial acceleration of shells 
in the outward direction. The energy of trapped radiation is thus 
spent to accelerate the outflow CM, and it is lost as a free en- 
ergy. A fraction (i*/i)~ 2 / 3 of the radiation energy survives till 
£* and contributes to the observed luminosity, 



to 



dt/frc 

At 



-2/3 



dt + U? Tee . 



The radiative efficiency (the ratio of to the explosion energy) 
is 

-2/3 

- ~ - - - * — I - ! [ ;in- + l] . (9) 



TcmMc 2 
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3 to 



Here, A = r rms0 /rcM is the initial amplitude of the fluctua- 
tions. 

In the case of Tqm > r„, we have t > an d all the radi- 
ated free energy will escape. Then, 

A 2 

(10) 



U free 



r CM Mc 2 2 ' 

If the initial Lorentz factor takes random values between 
r„ii n and r„i av (as assumed in most of the previous works), 



mm and T 
then 



A 2 = 



4-0 



3(1 + V) 2 ' 



(11) 



Since A 2 < 1/3 in equation (11), e does not exceeds 15%. 
If one further assumes that only 1/3 of the energy is actually 
radiated and the other 2/3 are stored as internal energy subject 
to adiabatic cooling, then the efficiency is reduced by a factor 
of 3 (2/3 of the dissipated energy is then spent to accelerate the 
CM). One thus arrives at <~ 5 % limit on e. 



C/fre 
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(6) 



M 



where M is the total mass of the outflow. This energy is avail- 
able for dissipation. From equation (5) one gets 



t/fr, 



^froo ( — 



-2/3 



(7) 



In Figure 1, we illustrate solution (5) by direct numerical simu- 
lations. 

Most of the free energy is dissipated at t ~ to- It should be 
compared with the time at which the outflow gets optically thin, 



Mar 



87TTO p C 2 r 2 



87TTO p C 4 r 3 



2 x lO^IT 



(8) 



where M is the mass outflow rate and L = Mc 2 T av is the ki- 
netic luminosity (assuming a spherical outflow). In the case of 

Tcm < T, w 180(r rms /r CM ) 1/5 i52 1/5 (Ao/3xl0 7 )- 1 /5, dis- 
sipation starts at to < t* . The free energy remaining in the out- 
flow by the transparency moment is [7 f * co = [7 f ? oc (i*/io)~ 2 / 3 , 
and most of the energy is dissipated at the optically thick stage. 
The radiation produced at t < t* is trapped in the plasma whose 
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FIG. 1. — Evolution of N = 3000 shells with initial separation 1 ms and 
r fluctuating according to equation (12) with A = 0.2 and To = 100. The 
shells begin to collide at to ~ 50 s. Solid curves show the average mass of the 
merging shells, m(f), the rms of their Lorentz factors, r rms (t), and the free 
energy of the system, C/f rce (t). Dashed lines display the self-similar solution 
(5,7). The deviations from the power-laws occur at late times when the number 
of shells is reduced and most of the free energy has been dissipated. 
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Numerical simulations of Kobayashi et al. (1997; see their 
Table 1) are in good agreement with equations (10,11), ex- 
cept for the case where tp — > and m, oc T^ 1 (their pa- 
rameter rj = — 1). One can show that in this special case 
rcM/r rms Rj v / V' ~~ * and e w l + y/tflni/) — > 1. The outflow 
dynamics is then non-linear (even though A = 1/V3 < 1), 
similar to the high amplitude case we study in § 3. 

3. HIGH- AMPLITUDE FLUCTUATIONS 

Suppose that T — 1 fluctuates with a log-normal distribution, 

r-i e-« 2 / 2 

In= ~ = M, P(0 = -7=— (12) 

1 - J- V 27T 

Here, ^4 measures the amplitude of the fluctuations. At A < 1, 
we have (T - r )/r w A£, so that T av = T and r rms /r av = 
A. At A > 1 we are in the high-amplitude regime. 

The possibility of high efficiencies at A > 1 can be under- 
stood when comparing two global characteristics of the out- 
flow: the CM Lorentz factor and the specific energy, n (in units 
of c 2 ). For definiteness, take a uniform initial distribution of 
mass. Then, 

CM "7r^' (13) 



ry = y P(0(r-1)<^. (14) 

It easy to show that at high amplitudes r\ 3> Tqm, i.e., the 
chaotic (free) component of the kinetic energy is much larger 
than the regular component. If the dissipation occurs in the op- 
tically thin regime and the emitted radiation is approximately 
isotropic in the CM-frame, then the momentum conservation 
law implies Tqm ~ const. When the dissipation is done, the 
final kinetic energy is about (Tcm — l)Mc 2 <C rjMc 2 , i.e., 
almost all the explosion energy has been radiated to infinity. 

In fact, in the high-amplitude case, the fluctuations start to 
dissipate very early, much before the transparency moment. 
The optically thick dissipation accelerates the CM, as discussed 
in § 2: the free energy gets transformed into the bulk motion of 
the outflow as a whole. The process is non-linear: the CM is 
substantially accelerated, in contrast to the low-amplitude case. 
It tends to reduce the initially big difference between Tqm and 
-q by increasing Tcm, while rj stays at the initial value, n = t]q. 

If at the transparency moment the ratio n/T cm is still high, 
then a high efficiency can be expected: the dissipation at the 
optically thin stage converts the difference between r\M c 2 and 
(Tom — l)Mc 2 into the observed radiation. Note also, that 
the radiation produced before the transparency moment is not 
completely lost. Even being cooled adiabatically, it contributes 
substantially to the outgoing luminosity. 

We now illustrate with numerical simulations. We generate 
a sequence of 3 x 10 3 thin shells with initial separation 1 ms. 
Their Lorentz factor fluctuates according to equation (12). The 
shells have equal mass m . Note that the results depend on 
the initial mass distribution (in contrast to the linear case), and 
rrii = mo is taken as a simple example. The duration of the 
central engine activity is 3 s. The results do not change sub- 
stantially if one takes longer activity: each 1 s portion of the 
outflow is causally disconnected from the other portions during 
the main emission time, t < 10 4 s. 



The transparency moment, t*, is roughly estimated by equa- 
tion (8) with r ~ Tcm- The transition from the optically thick 
to optically thin regime is treated in the simplest way: (1) If two 
shells merge at t < t», we assume that no radiation is emitted. 
Instead, the radiation is trapped and contributes to the kinetic 
energy of the new big shell. In other words, shell coalescence 
at t < t* proceeds with conservation of energy, n = const. A 
fraction (t*/t)~ 2 / 3 of radiation trapped at moment t survives 
till t* and contributes to the luminosity, L„. (2) At t > t*, 
only specific momentum conserves in the coalescence events. 
The energy released in the inelastic collision is radiated away 
isotropically in the rest frame of the newly formed shell. 

We consider models of two types. Model I: t* = 200 s and 
T in equation (12) is adjusted in such a way that Tcm ~ 100 
after the dissipation is finished. Model II: t* = 7 s and To is ad- 
justed to get the final Tcm ~ 300. In both models, the isotropic 
kinetic luminosity, L <~ 10 52 erg s _1 . 

The outflow forms by t = 3 s and after that it has a well 
defined CM with a velocity (3cm = Pc/E, where P is the to- 
tal momentum and E = (n + l)Mc 2 is the total energy of the 
outflow. (Note that the initial rj = rjo depends on the specific 
statistical realization I\ and it varies around the average ex- 
pected value given by eq. [14]). In Figure 2, we show examples 
of the non-linear evolution of rj and Tqm in Models I and II. 
At high A, the radiation is mostly produced by shells moving 
faster than the CM, and this results in the CM deceleration at 
t > i*. 

Then, we compute a sequence of models and find the effi- 
ciency, e = L^/ijoMc 2 , as a function of A (Fig. 3). At A < 1, 
we are in the linear regime of § 2. Shells begin to collide at 
to w yl _1 ro(Ao/c) and then evolve according to solution (5). 
Model II (Tcm ~ 300 > T*) is in perfect agreement with 
equation (10). The efficiency of Model I (r CM ~ 100 < T*) is 
reduced as a result of adiabatic losses during the optically thick 
stage. 
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FIG. 2. — Examples of fireball evolution with a strongly fluctuating Lorentz 
factor (A = 4 in eq. [12]) for Model I and II (see the text). The outflow is 
formed during the first 3 s. Dashed curves show the specific energy of the 
outflow, rj(t), and solid curves show the center-of-momentum Lorentz factor, 
Tcm- Dotted lines mark the transparency moments, t„, for the two models. 
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At A > 1, the efficiency increases up to ~ 83 % in Model I, 
and up to 96 % in Model II. The global parameters of the out- 
flow show substantial variations depending on the particular re- 
alization of the initial Tj. Figure 3 shows e averaged for many 
realizations and its standard deviation. 

Figure 4 shows examples of the generated light curves. Each 
coalescence event produces a pulse of a standard shape corre- 
sponding to a thin instantaneously radiating shell. The observed 
pulse has a width Ai Q b s ~ t/T 2 , where t is the time at which 
the radiation leaves the shell. Although the initial fluctuations 
are Gaussian, the non-linear dissipation generates a highly cor- 
related signal which is a mixture of quiescent periods and pe- 
riods of strong activity. Such a behavior is observed in GRBs 
and this special feature is naturally explained by the non-linear 
model. 

There is another special feature of the non-linear regime. The 
effective temperature of fluctuations in the CM-frame is rela- 
tivistic, and an additional emission mechanism appears: inverse 
Compton scattering (IC) on the bulk motions (see also Lazzati 
et al. 1999). This mechanism does not work at the optically 
thick stage since the radiation is trapped and follows the mo- 
tions of shells. However, when the outflow gets optically thin, 
the radiation can propagate and promote momentum exchange 
between the shells without direct collisions (like Silk damping 
in the early Universe). This effect should be accounted for in 
future. We expect that it will not change crucially the efficiency 
because: (1) Anyway, the free energy is radiated away, what- 
ever dissipation mechanism works. (2) Photon exchange does 
not crucially enhance the rate of dissipation because direct col- 
lisions also occur with relativistic velocities in the CM-frame. 

The IC by bulk motions can provide an observational test 
for the model. In particular, it can explain the high energy tail 
sometimes observed in GRB spectra. The photons emitted by 
internal shocks (with energy ~ 10 — 1000 keV) are boosted in 
energy by a factor of (T lms /T cm) 2 , resulting in GeV emission. 

The study of fluctuating fireballs in this paper was limited to 
the case of white fluctuations. Fluctuations with an arbitrary 
spectrum will be studied in a future paper (Beloborodov 2000). 
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FIG. 3. — Efficiency, e, as a function of the fluctuation amplitude, A. Open 
circles show the results in Model I (with final Tcm ~ 100) and filled circles 
show Model II (with final Tcm ~ 300). Bars show the standard deviation 
(where it is larger than the symbol size). Line displays e = A 2 /2 expected in 
the optically thin linear regime (see eq. [10]). 
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FIG. 4. — Example light curves in Model I with A = 1 and A = 4. 
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